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Abstract 

We investigate nonlocal coherent transport in ballistic four-terminal 
Josephson structures (where bulk superconductors (terminals) are con- 
nected through a clean normal layer, e.g., a two-dimensional electron gas). 

Coherent anisotropic superposition of macroscopic wave functions of 
the superconductors in the normal region produces phase slip lines (2D 
analogs to phase slip centres) and time-reversal symmetry breaking 2D 
vortex states in it, as well as such effects as phase dragging and magnetic 
flux transfer. The tunneling density of local Andreev states in the normal 
layer was shown to contain peaks at the positions controlled by the phase 
differences between the terminals. 

We have obtained general dependence of these effects on the control- 
ling sup ercurrent /phase differences between the terminals of the ballistic 
mesoscopic four-terminal SQUID. 
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1 Introduction. 



The multiterminal Josephson junctions [jl], |[ generalizes the usual (two-terminal) 
Josephson junctions |3| to the case of weak coupling between several massive su- 
perconducting banks (terminals). Compared with two-terminal junctions, such 
systems have additional degrees of freedom and the corresponding set of control 
parameters, preset transport currents and (or) applied magnetic fluxes. As a 
result, the current- or voltage-biased and the magnetic flux-driven regimes can 
be combined in one multiterminal microstructure. 

One of the realizations of multiterminal coupling is a system of short dirty 
microbridges going from a common center to separate massive superconduc- 
tors. The theory of this kind of multi-terminals was derived in |B| within the 
phenomenological Ginzburg-Landau scheme ( Aslamazov and Larkin model B ) . 
This approach is valid for temperatures T near the critical temperature T c and 
for the local case when the characteristic spatial scale is larger then the coher- 
ence length £o ~ %Vf/T c . The stationary states and the dynamical behaviour of 
the microbridge type multitcrminals were studied for different microstructures, 
four-terminal SQUID controlled by the transport current, weakly coupled su- 
perconducting rings (see review of theoretical and experimental results in R W) . 

The Josephson effect in mesoscopic weak links with direct conductivity (S- 
N-S junctions , ballistic point contacts) exhibits specific features 0,pQ] which 
are absent in conventional dirty microconstrictions near T c [o. As in normal 
metal mesoscopic structures JTT| , the electrodynamics of supercurrents in meso- 
scopic Josephson junctions is nonlocal. Supercurrent density depends on the 
spatial distribution of the superconducting order parameter in all points of the 
mesoscopic weak link region. The coherent current flow is carried by the An- 
dreev states formed inside the weak link. Nonlocal nature of mesoscopic 
supercurrents was demonstrated by Heida et al.p^, investigating the meso- 
scopic S-2DEG-S (superconductor-two dimensional electron gas - superconduc- 
tor) Josephson junctions. They measured 2$o periodicity of the critical current 
instead of the standard <E>o (^o = hc/2e is the magnetic flux quantum). Theory 
of this effect was developed in Refs.jlJ, 

The current level of nanofabrication technology made it possible to realize 
multi-terminal mesoscopic Josephson junction similar to the 2-terminal junction 
studied in |13| , The microscopic theory of the mesoscopic ballistic Josephson 
multi-terminals was derived in Rcf . [l6| . It is valid for arbitrary temperatures 
< T < T c and describes the nonlocal coherent current states in the system. 
The effects of nonlocal coupling, such as phase dragging and magnetic flux 
transfer were obtained in Ref. |L7| . 

In the present paper we continue study of quantum interference effects in 
mesoscopic multitcrminals, which are related to the nonlocality of weak cou- 
pling. The paper consists of two parts. In first part (Article I) the effects of 
nonlocal coupling in mesoscopic multiterminal structures are studied. The gen- 
eral properties of Josephson multiterminals are described in Section 2. Section 
3 gives the results concerning the current distribution and density of states in- 
side weak link. In Section 3 we study specific for mesoscopic case properties of 
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Figure 1: Mesoscopic four-terminal Josephson junction with "parallel" imple- 
mentation of the supercurrents. The four bulk superconducting regions, S1...S4, 
are weakly coupled through the thin layer of normal metal (2DEG), represented 
by the shaded area. 



four-terminal SQUID. In the second part (Article II) a superconducting phase 
qubit based on mesoscopic multiterminal junction is proposed and investigated. 



2 Mesoscopic Four- Terminal Junction. 

2.1 System description 

In a mesoscopic 4-terminal junction, the bulk superconductors (terminals) are 
weakly coupled to each other through a clean two-dimensional normal metal 
layer (2D electron gas) as it is shown in Fig.l. The pairs of terminals can be 
incorporated in bulk superconducting rings or in circuits with preset transport 
currents. In Fig. 2 we show two such configurations. The first one (Fig. 2a) 
presents two superconducting rings, each interrupted by a Josephson junction, 
which are at the same time weakly coupled to each other. The second con- 
figuration (Fig. 2b), combines a current (or voltage) biased junction and a flux 
driven junction in the ring. We call this configuration the 4-terminal SQUID 
controlled by the transport current. 

The state of the i-th terminal S, (i=1...4) is determined by the phase ipi of 
the complex off-diagonal potential Aoexp(i^). Superconducting banks induce 
the order parameter 4' in the normal metal region (shaded area in Fig.l). Inside 
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Figure 2: Superconducting microstructures based on mesoscopic four-terminal 
Josephson junctions, (a) Two weakly coupled superconducting rings, (b) Meso- 
scopic four-terminal SQUID. 



this mesoscopic, fully phase coherent weak link, the supercurrent density j(p) 
at point p nonlocally depends on the values of the induced order parameter 'J 
at all points p '. In its turn, the order parameter ^f(p) depends on the phases 
ipi. The total current Ii, flowing into i-th terminal depends on the phases tpj of 
all the banks and has the form [flq] : 



7*"A . / <fi 



tanh 



A cos(^) 



2T 



(1) 



In the case of two terminals Eq. (1) reduces to the formula for ballistic point 
contact [[l0| with 712 equals to Sharvin's conductance. 

Expression (1) corresponds to the case of small junction, when linear di- 
mensions of the N-layer are smaller than the coherence length £ ~ Tivp/A.^ 
(for the case of arbitrary junction's dimensions see Ref . Jl6| ) . We are focusing 
here on the small junction case because the effects of nonlocality are most pro- 
nounced in this situation. The geometry dependent coefficients 7^ denote the 
coupling between the partial Josephson currents in ballistic two-terminal Si-S^ 
weak links. 

Eq. (1) is simplified when T ■ 



0, or when T ss T c . In T = limit, it becomes 



^Aq(O) 



Sign 



<fi - ifj 



Near T c on the other hand, the order parameter is small, Ao 
write 

\2 4 



U, = — r-= — 2^ sm{<pi 



4eT c 



(2) 

and one can 

(3) 



Equations (2) and (3) are qualitatively similar, differing by the magnitude of 
critical currents and by the shape of current-phase dependencies (sin(^)sign[cos(£)] 
and sm(<^)). For definiteness, in the following we will consider the case of T ~ T c , 
keeping in mind that the results hold qualitatively at low temperatures as well. 
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For the Josephson coupling energy of the junction Ej , which relates to the 
supercurrents U (3) through Jj = (2e/K)dEj /d<fi, we have 



h 7rAo(T) 2 
2e 4eT c 



5^7ife[l ~ cos (<pj -<pk)]. 



(4) 



Expression (1) for supercurrents Ii looks similar to Buttiker's multiprobe 



which relates the currents to the voltage drops between terminals in mesoscopic 
normal metal multiterminal system. The similarity reflects the above mentioned 
nonlocality of mesoscopic transport on the scale of £t ~ hvp/T (in the ballistic 
limit we are considering). The essential difference between (1) and (5) is that 
unlike Josephson currents of (1), the normal currents of (5) can flow only out of 
equilibrium; while the current-phase dependence in (1) is periodic, the current- 
bias dependence of (5) is certainly not. 

2.2 Circuit implementations of 4-terminal junction. Non- 
local weak coupling. 

The current-phase relations (3) determine the behaviour of the system in the 
presence of the transport currents and/or the diamagnetic currents induced by 
the magnetic fluxes through the closed superconducting rings. It is necessary 
to distinguish two types of circuit implementation of the mesoscopic 4-terminal 
junction 0. The first one, is the "crossed" or "transverse" implementation, 
when the total current in one circuit goes in and out through one pair of oppo- 
site banks in Fig.l and in the second circuit - through the other pair. In the 
"parallel" implementation, shown in Fig.l, the currents I and J flow through 
the pairs of adjacent banks. In this case, nonlocal coupling of currents inside the 
mesoscopic N-layer results in peculiar effect of "dragging" of the phase differ- 
ence between one pair of terminals by the phase difference between another pair 
of terminals P7J . In the following, we consider the "parallel" implementation 
and study the manifestations of the phase dragging effect. 

The coefficients in (3,4) depend on the geometry of the weak link (the 
shape of the N-layer) and on the transparency of S-N interfaces. In general we 
have jij — jji and = 0. For the case of parallel implementation, the elements 
712 and 734 are related to the critical currents of the individual sub-junctions 
S1-S2 and S3-S4 respectively. The matrix 



describes the coupling between these two junctions. Properties of the system (in 
particular the existence of the phase dragging) qualitatively depends on whether 
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(5) 
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det(% oup i) equals to zero or not (see Appendix). In case of a conventional non- 
mesoscopic 4-terminal Josephson junction the coefficients 7^ factorize, 7^ ~ 
(1/Ri)(l/Rj), where Ri are the normal resistances of dirty microbridges 0. 
This yields det(j coup \) = 0, which we call local coupling. On the other hand, 
in a mesoscopic system, even in a completely symmetric case of a x a square 
N-layer and ideal transparency (D = 1) of N-S$ interfaces, the coefficients 7.^ 
are given by ]l7| : 

(y/2 1 \ e 2 p F a ( 1 \ 

712 = 734 = 70, 7cou P i = 7o^ 1 70 = -^-^^-^, (7) 

with det(% oup i) 7^ 0. In more general case than the completely symmetric one 
(Eq. (7)), we can write 7^ in the form 

734 = K712 , 7cou P i = ' ( 8 ) 

This corresponds to a square N-layer, with different transparencies for junctions 
S1-S2 and S3-S4 and/or different width of the superconductor banks connected 
to the normal layer. In our numerical calculations wc will use the simple form 
(7), i.e. K=l,p=V2,q=l. 

2.3 Current-phase relations. The phase dragging effect. 

Let us introduce new variables: 

<P2 ~ Pi = , <f3 — ipi = (p 

a~[3 = x , a + /3 = 7 . (9) 

Without loss of generality, we can choose the phase 7 equal to zero . <f>j = 0). 
For the circuit implementation shown in Fig.l, we have 

I = h = -h, J = h = -h. (10) 
In terms of phase differences (9) the currents / and J have the form 



I = suit 



J = k sin i 



+ q) sin - cos - + (p - q) cos - sin - 



cosx, (11) 



(p + q) sin - cos - + (p - q) cos - sin - 



cos X . (12) 



All 7's (8) are normalized by 712 and the currents I, J are measured in units of 
I =7T7 12 A (T) 2 /4eT c . 

From the current conservation (Eq.(10)), it follows that the phase \ m Eqs. 
(11) and (12) can take only two values, or tt. Minimization of Ej (4) with 
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Figure 3: Phase diagram for the phase difference \ m the (9, ^)-plane. Solid 
line separates the regions with \ — and x = n - The dashed region is absent 
in the case of local coupling. 



respect to \ also gives x = or tt, depending on the equilibrium values of 9 and 
<j> (see Appendix): 



cosx = sign 



(p + q) cos - cos - 



(P- 



q) sin - sin - 



(13) 



The current-phase relations (11) and (12) with the condition (13) are invari- 
ant under the transformation 9 — > 9+2irn, and <f> — > tfi+2irk. The 2ir periodicity 
of observable quantities is sustained by the "hidden" variable phase x- I n Fig. 3 
the phase diagram for \ m the (9, 4>) plane is presented. The solid line separates 
the regions with % = and \ — n - When the state of the system (9, <fr) crosses 
this line, a jump in x occurs. Corresponding jumps take place in current-phase 
relations (11) and (12). The current 1(9) (11) is shown in Fig. 4 for several 
values of the phase <f>. Note that the function 1(9) has jumps, which for 0^0, 
are located not at 6 = ±7r, as they would be in conventional 4-terminal junc- 
tions. The jump in x means the slippage of the phase 9 (or <fi). In the case of 
two-terminal or conventional 4-terminal junction the phase-slip events occur at 
phase difference equal to n(2n + 1), n — 0, ±1, ±2, ... In one-dimensional struc- 
tures slippage of the phase occurs at phase-slip centers (PSC), i.e. points where 
the order parameter equals to zero. In our case of 2D mesoscopic 4-terminal 
weak link, the analog of the PSC are phase-slip lines in normal metal region. 
They appear when the state of system (9, <f>) belongs to the dashed region in 
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Figure 4: Current-phase relations 1(8) for different values of <fi. 

Fig. 3. This region, which is absent in the local coupling case (it actually coin- 
cides with lines 8 — ir, <fi = ir), we call "frustrated" region for phases 8 and (j>. 
For states inside this region, the distribution of the supercurrent in the weak 
link contains 2D vortex states (see below). 

Nonlocal weak coupling leads to the phase dragging effect [jl7|. One notices 
that if p q then putting 6 = in (11) results in a nonzero value of the current 
/ 



This current is absent in conventional 4-terminal junctions or mesoscopic four 
terminal junctions with crossed implementation at which p = q (i.e. cfei(7 coup i) = 



Similarly if we set / = Oin (11), we find a nonzero solution for 8, which again 
vanishes whenp = q. This solution (= 8d) is a function of <f> and is plotted versus 
4> in Fig. 5. The influence of the phase of one side of the mesoscopic 4-terminal 
junction on the phase of the other side is what we call phase dragging effect. This 
effect is one of important characteristics of the with parallel implementation. 

In general current-phase relations are asymmetric, I(—0) ^ —1(8), unlike in 
conventional cases. In another words, the presence of phase difference <fi on the 
terminals S4-S3 breaks the time reversal symmetry for Josephson junction S1-S2. 
It also follows from expression (11) that 1(8) is not only a function of \(f>\, as in 
conventional junctions, but also depends on the sign of 0. The phase dragging 




(14) 



0). 




Figure 5: The dragged phase Od between terminals S1-S2, at zero transport 
current /, as a function of the phase difference cf> between the other pair of 
terminals S3-S4. 



has the analogy in the normal metal mesoscopic multiterminals, described by 
formula (5); the normal current flowing through one pair of terminals induces a 
voltage difference between the other ones Q . 

3 Current distribution and local density of states 
inside the mesoscopic weak link 

The coupling through the normal layer determines the behaviour of the Joseph- 
son weak links S1-S2 and S4-S3. On the other hand, the properties of the 
normal layer itself depend on the phase differences 8 and <fi on the junctions. 
The phases 9 and <fi can be controlled by external magnetic fluxes through the 
rings (Fig. 2a). In this section we present the results of numerical calculations 
for current density distribution j(p) and density of local Andrccv levels N(e) in- 
side the mesoscopic 4-terminal weak link. The expressions for j(p) and N(e) as 
functionals of {(fi, 1^2,^3, f^} were obtained in Ref. [|l6| by solving Eilenberger 
equations |l9| ]. 

Fig. 6 illustrates the effect of phase dragging. Two sets of phases (6 = —0.42, 
4> = 7r) and (9 = 0.42, (j> = ir) correspond to zero value of the current / (II) 
(see Fig. 4) and opposite directions of the current J (12). In the absence of the 
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Figure 6: Distribution of the current density inside the normal layer for phase 
<f> = 7r and two values of the phase 9 at which the current I = (Fig. 4). (a) 
9 = -0.42, (b) 9 = 0.42. 
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Figure 7: The vortex-like distributions of the current inside the weak link when 
and </) are inside the frustrated region: (a) 6 = <fi = n — 0.2, (b) = <f> = 7T + 0.2 



current from terminal Si to terminal S2, the phase difference on the junction 
S1-S2 exists. 

When the phases 9 and <j> lie in the "frustrated" region of the diagram Fig. 3 
(dashed area), the current distribution j(p) contains 2D vortex states. They are 
shown in Fig.7 for states (6 = tt-0.2, = tt-0.2) and (6 = tt + 0.2, <j> = tt + 0.2). 
In both cases, the order parameter ^(p) vanishes along the diagonal x — y and 
its phase drops by tt when crossing this 2D phase-slip line. 

The Andreev scattering processes on the Sj-N interfaces lead to the appear- 
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Figure 8: Density of states, iV(e), averaged over the normal region for different 
values of 9 and <j>. 



ance of the energy levels with energies e inside the gap A , |e| < A , in the 
normal metal. The local density of electron states in the normal layer is given 
by the formula 

Af(e, p) = N(0) < Re g(cu = -it, p, v F ) > VF . (15) 

( g(oj, p, vp)) is Eilenberger Green's function). We have studied the dependence 
of the density of states, averaged over area of the N-layer, N(e), on the phases 9 
and <fi. This tunneling density of states can be measured by scanning tunneling 
microscope. It contains the spikes with intensity and position on the energy 
axes controlled by the phases 6 and 4>. The results are shown in Fig. 8 ( the 
(5-function singularities in N(e) are smeared by introducing a small damping 
r = O.OIAq). 
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4 Mesoscopic Four- Terminal SQUID 



In this section we consider the four-terminal SQUID configuration (Fig. 2b). 
Conventional 4-terminal SQUID has been studied in detail in Ref. ||, wherein 
the steady states domain and dynamical properties of the system were calcu- 
lated. Here we are interested in the specific features of the mesoscopic case 
reflected in the current-phase relations (11,12). As we have seen in the previous 
section, the nonlocal coupling (p ^ q) leads to the phase dragging effect. This 
dragged phase can induce a transferred magnetic flux in the ring which depends 
on the transport current. On the other side, the magnetic flux state in the ring 
influences the behaviour of the Josephson junction in the current circuit. 

When the terminals 3 and 4 are short circuited by a superconducting ring 
with self- inductance L, the phase 4> is related to the observable quantity, mag- 
netic flux threading the ring $, <fi — ^ $. The current J circulating in the ring 
is given by J = (<& e — $)/L, where $ e is the external magnetic flux threading 
the ring. From (11) and (12) we have 



/ = sin ( 



9 $ 9 $ 

(p + q) sin - cos - + (p - q) cos - sin — 



cosx, 



(16) 



$ e - $ 1 

= sin $ + - 

C K 



$ $6 

(p + q) sin - cos - + (p - q) cos - sin - 



cosx, 



(17) 



where fluxes $,$ e are measured in units h/2e, C = (2e/h)LIoK is the dimen- 
sionless self- inductance. The parameter k = 734/712 is the ratio of the critical 
currents of the sub-junctions 3 — 4 and 1 — 2. The limiting cases of n — > 00 and 
k — ► correspond to the autonomous SQUID and the current biased Josephson 
junction, respectively. 

The transport current I and the external flux $ e are the external controlling 
parameters. The corresponding Gibbs potential for the 4-terminal SQUID takes 
the form 

U(<$>, 9; I, $ e ) = V — ' -19- cos(6>) - k cos($) 

/ $ 9 $\ 

-2 I (p + q) cos - cos - - (p - q) sin - sin - I cos(x) (18) 

The last three terms in Eq.(18) are the Josephson coupling energy (4) in terms 
of variables 9 and x- The minimization of U with respect to x gives the 
expression (13) for cos(x), with (f) replaced by <&. At given values of the control 
parameters / and $ e , the relations (16) and (17) (together with Eq.(13)) de- 
termine the set of possible states of the system {9, <E>}, among which we should 
choose those that correspond to the local minima of the potential U, Eq.(18). 
Let us consider the effect of the magnetic flux state of the ring on the behaviour 
of the current driven junction. The critical current of the junction, I c , depends 
on the applied magnetic flux <£> e . In the simplest case of small self-inductance 
C -C 1, we can neglect the difference between $ and $ e in expression (16). 
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Figure 9: The steady state domain for mesoscopic four-terminal SQUID in plane 
(I, $ e ) of the control parameters (solid line). Dashed line corresponds to the 
conventional four-terminal SQUIDs. 



The maximal value of the supercurrent / (16) (with $ replaced by $ e ) as a 
function of $ e , I ma x{& e ), is shown in Fig. 9. This curve determines the boundary 
of the steady states domain in the (I, $ e ) plane. The function I m ax(& e ) is 
2tt periodic, but due to the terms proportional to p — q in Eq.(16), it is not 
invariant under the transformation $ e — ► — <£> e . The symmetry is restored if we 
simultaneously change <& e on — <!> e and I on —I. Note, that in conventional case 
(p = q) the boundary of the steady state domain I ma x(^ e ) is symmetric with 
respect to the axes (J, $ e ) (dashed line in Fig. 9). Thus, the critical current I c 
in the transport current circuit, for a given direction of the current, depends on 
the sign of the magnetic flux in the ring. For finite values of self-inductance C 
equations (16) and (17) must be treated self-consistently. The critical current 
I c as function of C is shown in Fig. 10 for two values of external flux, <I> e = 
and $ e = 7T. 

Outside the steady state domain, the stationary solutions for (8, $) are ab- 
sent and system goes to the nonstationary resistive regime. The simple gener- 
alization of Eqs. (16,17) in the frame of the heavily damped resistively shunted 
junction (RSJ) model || leads to equations (see 0): 



d6 
~dt 



I - sin e 



9 $ 8 $ 

(p + q) sin - cos — + (p - q) cos - sin — 



cosx, 



(19) 



~dt 



-sin<i>- 



$ 8 $6* 
(p + q) sin — cos - + (p - q) cos — sin - 



cosx, (20) 



13 



0.5 

L 



Figure 10: The critical current, I Cl between the superconductors Si and S2, as 
a function of C for <f> e = and ir. 



dx 
dt 



(p + q) cos - cos - 



sin(x) 

They can also be presented in a form 

$ = 



[P - q) sin - sin - 



dU 

de 



d<S> 



X 



ldU 

"2 ax' 



(21) 



(22) 



where potential U is defined in Eq.(18). The voltages between different terminals 
are related to the time derivatives of the phase differences 



Vn = e, ^ 34 -<i), l(v 13 + v 24 ) = x- 



(23) 



The time and the voltage are measured in the units of e/Jo and hIo/2e 2 re- 
spectively. Note that, inspite of the equilibrium state, the dynamical variable 
X relates to an observable quantity. It's time derivative determines the voltage 
between the ring and the transport circuit. The features of the dynamical be- 
haviour of the mesoscopic 4-terminal SQUID are again affected by the terms 
proportional to (p — q), i.e. by nonlocal coupling. The current- voltage charac- 
teristics in the transport channel, V(I), (the time averaged voltage V21 (23)), 
can be obtained by the numerical solution of the coupled system of nonlinear 
differential equations (20-22). As well as critical current I c , the voltage V(I) in 
applied magnetic flux $ e depends on the sign of the <E> e , i.e on the direction of 
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Figure 1 1 : The flux induced inside the ring as a function of the transport current 

I. £= l,$ e = 

the external magnetic field. Full dynamical description of the mesoscopic four- 
terminal SQUID requires more rigorous approach than RSJ model, and will be 
the the subject of separate investigation. In accordance with stationary (16,17) 
or dynamical equations (19-21) for and $ the opposite effect for influence of 
transport current circuit on the flux states in the ring takes place. In particular, 
the current I produces the flux $ in the ring even in stationary case and in the 
absence of external flux $ e . This effect is proportional to (p — q) and is absent 
in the conventional case. In Fig. 11 we plot the dependence of the induced in 
the ring magnetic flux <E> on the transport current I in the case $ 6 = 0. 

Special interest presents the existence of the bistable states in the system 
described by the potential (18). We emphasize, that in contrast to the usual 
SQUID, bistable states occur for any inductance C, even for C < 1 ||. We will 
analyze the dependence of these states on the control parameters I and $ e in 
article II, when designing of the four-terminal qubit will be studied. 

5 Conclusions 

We have demonstrated that in ballistic four-terminal Josephson junctions co- 
herent anisotropic superposition of macroscopic wave functions of the supercon- 
ductors in the normal region produces formation of phase slip lines (2D analogs 
to phase slip centres) and time-reversal symmetry breaking 2D vortex states 
in it, as well as such effects as phase dragging and magnetic flux transfer. We 
have calculated the phase-dependent tunneling density of Andreev states in this 
region as well. 
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A degree to which the nonlocality of mesoscopic transport is manifested, de- 
pends on the characteristics of the system and is most pronounced in the ballistic 
case [ p0| . Ballistic four-terminal junctions considered here demonstrate several 
specific effects absent in the diffusive limit || ||, 0]: the phase dragging, time- 
reversal symmetry breaking (1(6) ^ I(—8), Eq.(14)), and the vortex formation. 
The latter can mimic the behaviour of SNS junctions with unconventional su- 
perconductors pl| . It has indeed the same origin in direction-dependent phase 
of the superconducting order parameter induced in the normal part of the sys- 
tem, though not due to the intrinsic phase difference between different directions 
in a superconductor. This actually allows us more freedom in controlling the 
behaviour of the junction, which will be exploited in the qubit design based on 
such a junction in the following paper. The time-reversal symmetry breaking 
can be also used for direction-sensitive detection of weak magnetic fluxes. 

It will be instructive to investigate the role played by finite elastic scattering 
in the system and look for the analogs of zero bound states, found at sur- 
faces/interfaces of unconventional superconductors (for a review see [p2| ). This, 
as well as vortex dynamics in the system, will be the subject of our further 
research. 
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6 Appendix. Junction with arbitrary 7's 

The Josephson energy of the mesoscopic Four Terminal Junction, normalized to 
(h/2e)(irA (T) 2 /4eT c ), is expressed by 

Ej = -712 cos - 734 cos 4> + E coup i (24) 

with the coupling energy -E COU pi given by 



-Ecoupi = -713 cos h x -714 cos 



2 



"723 cos ( — h x ) - 724 cos [ — h x 



= — (A cos x + B sin x) (25) 

where 



A = (713 + 724) cos ( - ) + (714 + 723) cos 
B = (713 - 724) sin ( ^r 1 - ) + (714 - 723) sin ( "— — ) (26) 
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Minimizing with respect to x, we find the minimum to be 



Scoupl = -V^+B 2 , 

After some manipulations we find 

-^coupl — 



COS 



.4 



VA 2 + B 2 



7l3 + 723 + 7l4 + 724 



(27) 



(28) 



+2(7i37i4 + 723724) cos0 + 2(713723 + 714724) cos6» 

+2(713724 + 714723) cos 8 cos 4> - 2(7x3724 - 714723) sin6>sin</> ] 1/2 



The last term in the bracket in (|2£|) vanishes when det(^ coup i) = 0. In that case 
the current 1(8, <f>) will be zero at 8 = 0. On the other hand, if det(j coup i) 7^ 
then 1(8, <fi) ^ when 8 — 0. This is a signature of the phase dragging effect. 

In a four terminal junction with micro-bridges near T c one has 7^ ~ 1/RiRj. 
In that case the last term in (|2q ) will vanish and -E C oupi factorizes 



coupl 



1 



1 

S2" 



4cos 2 (0/2) 
R1R2 



1 



1 
i?4 



4cos 2 (<^/2) 



R3R4 



In particular when R\ = R2 and i?3 = i?4 we find 

4 



E, 



coupl 



R1R3 



COS ■ 



COS ■ 



(29) 
(30) 



which is what one obtains from Ginzburg-Landau calculation. 

In a mesoscopic four terminal junction with parallel implementation on the 
other hand, we have 713 = 724 and 714 = 723. This leads to B = and therefore 



-Ecoupi = - \A\ , cos x = sign(A) 



(31) 



which gives \ — or tt. Notice that in the general case of (|27j), \ can take other 
values than and %. 
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